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Abstract. Evidence for a large scale flow of low density 
gas onto the Cepheus A young stellar cluster is presented. 
Observations of K-band near-infrared and multi-transition 
CS and N2H + millimeter line emission are shown in re- 
lation to a sub-millimeter map of the cool dust around 
the most embedded stars. The near-infrared emission is 
offset from the dust peak suggesting a shift in the loca- 
tion of star formation over the history of the core. The CS 
emission is concentrated toward the core center but N2H + 
peaks in two main cores offset from the center, opposite 
to the chemistry observed in low mass cores. A starless 
core with strong CS but weak N2H" 1 " emission is found to- 
ward the western edge of the region. The average CS(2-1) 
spectrum over the cluster forming core is asymmetrically 
self-absorbed suggesting infall. We analyze the large scale 
dynamics by applying a one-dimensional radiative trans- 
fer code to a model spherical core with constant temper- 
ature and linewidth, and a density profile measured from 
an archival 850 /im map of the region. The best fit model 
that matches the three CS profiles requires a low CS abun- 
dance in the core and an outer, infalling envelope with a 
low density and undepleted CS abundance. The integrated 
intensities of the two N2H + lines is well matched with a 
constant N2H + abundance. The envelope infall velocity is 
tightly constrained by the CS(2-1) asymmetry and is sub- 
sonic but the size of the infalling region is poorly deter- 
mined. The picture of a high density center with depleted 
CS slowly accreting a low density outer envelope with nor- 
mal CS abundance suggests that core growth occurs at 
least partially by the dissipation of turbulent support on 
large scales. 
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kinematics and dynamics - ISM: molecules - ISM: abun- 
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1. Introduction 

Most stars, particularly massive stars, form in groups (e.g. 
|Carpenter| 2000). It is therefore essential to study clus- 
ter forming regions in order to understand more com- 
pletely the way in which the majority of stars are formed. 
Isolated low mass star formation occurs via the nearly 
isothermal free-fall collapse of a dense molecular cloud 
core, followed by the evolutionary phases Class 0, I, II 
and III objects (e.g. lEvansl 19991. However, the applica- 
bility of this paradigm to the formation of massive stars 
is debated ( |Caray fc Lizano| 1999): for example, massive 
stars begin burning hydrogen and reach the main sequence 
while still accreting matter from the surrounding proto- 
stellar envelope and they can also develop strong winds, 
both of which will strongly affect the physical conditions, 
structure and chemistry of their surroundings. Due to the 
shape of the IMF and the fact that they evolve faster, 
massive protostars are rarer (and therefore more distant 
on average) than low mass protostars. Consequently fewer 
Class massive protostar counterparts have been studied 
in detail. It is only recently that catalogs of high-mass 
protostellar objects have been made fe.g. lSridharan et al.l 
2002). 

The molecular cloud core Cepheus A East (hereafter 
Cep A) is a nearby site of massive star formation QSargent| 
1977) located in the Cepheus OB association at a dis- 
tance of 725 pc (jBlauuw et al.l 1959 ). The far-IR lumi- 
nosity is 2.4xl0 4 Lq (Eva ns~et al.l 1981), correspond- 
ing to a small cluster of B stars. Cep A harbors one 
of the first molecular bipolar outflow sources discovered 
dRodriguez et al.| 1980) . Higher spatial resolution CO ob- 
servations showed the outflow to be extremely complex, 
and it was termed quadrupolar UTorrelles et al. 1 1993). The 
fastest components of this outflow are bipolar and oriented 
northwest-southeast jRodrigu ez et al.| 1980), perpendicu- 
lar to the low velocity CO structure. The slower and more 
extended component has been interpreted as the divert- 
ing and redirecting of the main outflow by the interaction 
with interstellar high-density gas, seen in NH3 lines by 
ITorrelles et al.l f!993). Ultra compact H11 regions and a 
diffuse thermal dust emission source have been identified 
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from 20 /jm maps and 6 cm low-resolution VLA observa- 
tions by Bcichma n et al.l (1979). Seven ionized hydrogen 
complexes lie in "strings" that form a "Y" tilted to the 
east ( |Hughes &: Wouterloot| 1984), the bifurcation point of 
which is coincident with the exciting source of the molecu- 
lar outflow. A cluster of compact radio sources have been 
identified as pre-main-sequence stars by Hughes (1988) 
due to their variability and the presence of OH and H2O 
maser emission. From subsequent ammonia VLA obser- 
vations, ^H^E^OiLl (1993) proposed that one of these 
radio sources, HW2, is a ~10-20 M protostar. The larger 
core surrounding the cluster has temperature 35 K and 
mass 200 - 300M Q |Moriarty-Schieven et aT] (1991). On 
the basis of its protostellar content, high luminosity and 
low temperature, and following the bolometric tempera- 
ture definition of lChen et al.l (1995L the Cep A core may 
be considered a high mass Class source. 

In order to examine the properties of this young, 
massive cluster forming region, we obtained near-infrared 
and millimeter wavelength multi-transition CS and N 2 H + 
data. The observations are discussed in Sj2 We derive the 
density profile from 850 /im continuum measurements and 
fit core-averaged spectra using a 1-D radiative transfer 
model in fJ3 Our results indicate large CS depletion in 
the central core and an outer undepleted accreting layer. 
We discuss these results and conclude in 21 

2. Observations 

2.1. Millimeter Data 

Observations were made with the 10 antenna Berkeley- 
Illinois-Maryland array 1 (BIMA) for two 8 hour tracks 
in CS(2-1) in April and May 1998 and one 8 hour track 
in N 2 H+(l-0) in May 1998, all in C-array. A seven field 
hexagonal mosaic was made with phase center, a(2000) = 
22 h 56 m 18 i !9, <5(2000) = 62°01'42'.'6. Amplitude and phase 
were calibrated using 5 minute observations of 2322+509 
interleaved with each 25 minute integration on source. The 
calibrator flux was 0.66 ± 0.22 Jy based on observations 
of Uranus during the middle of each track. The correlator 
was configured with two sets of 256 channels at a band- 
width of 12.5 MHz (0.15 km s -1 per channel) in each side- 
band and a total continuum bandwidth of 800 MHz. Data 
reduction was carried out using standard procedures in 
the MIRIAD package. The final maps covered a hexago- 
nal region ~ 3.'5 x 3'5 region at ~ 9" x 7" resolution. 

Complementary single-dish maps of the same lines 
were made at the Five College Radio Astronomy Obser- 
vatory 2 (FCRAO) 14 m telescope in December 1999 using 

1 The BIMA array is operated with support from the Na- 
tional Science Foundation under grants AST-9981308 to UC 
Berkeley, AST-9981363 to U. Illinois, and AST-9981289 to U. 
Maryland. 

2 FCRAO is supported in part by the National Science Foun- 
dation under grant AST-0100793 and is operated with permis- 



the SEQUOIA 16 beam array receiver and the FAAS back- 
end consisting of 15 autocorrelation spectrometers with 
1024 channels set to an effective resolution of 24 kHz 
(0.06 km s _1 ). The CS and N 2 H + lines were observed si- 
multaneously in frequency switching mode. The pointing 
and focus were checked every three hours on nearby SiO 
maser sources. Third order baselines were removed from 
the data and spectra coadded using the CLASS package. 
The resolution of the data is 50" and the final maps were 
Nyquist sampled over 8' x 8' centered on the BIMA phase 
center. A single spectrum of C 34 S(2-1) was also taken to- 
ward the map center using the same setup. 

The FCRAO data were combined with the BIMA data 
using maximum entropy deconvolution (using a gain for 
FCRAO at these frequencies of 43.7 Jy K" 1 ). The result- 
ing maps show the large scale structure observed in the 
single-dish map at the ~ 10" resolution of the interferom- 
eter map. These maps have an rms noise of 0.2 and 0.4 K 
per 0.5 km s" 1 channel, for CS (2-1) and N 2 H+ (1-0) re- 
spectively. All the flux is recovered in the combined map 
but features at intermediate scales, in the range 70"— 100", 
may be poorly represented (e.g. William s et al.l 2003 h 

Observations of CS(5-4), CS(7-6) and N 2 H+(3-2) 
were made at the Heinrich Hertz Telescope 3 (HHT) in 
November 1999. The data were taken using the SIS-230 
and SIS-345 receivers and AOS backend (2048 channels, 
48 kHz resolution) in on-the-fly (OTF) mode. Pointing 
and focus were checked using observations of Saturn and 
Orion IRc2. The final maps, made from several OTF maps 
in orthogonal scan directions, were coadded and first order 
baselines removed using the CLASS package. The resolu- 
tion of the observations is 31", 22", and 27" for CS(5-4), 
CS(7-6) and N 2 H+(3-2) respectively and the final maps 
covered 3f5x3.'5, with an rms noise of 0.4 K per 0.5 km s _1 
channel. 

To compare the near-infrared and millimeter wave- 
length line data with the cool dust emission around the 
most embedded (Class counterpart) stars in the clus- 
ter, we downloaded archival SCUBA observations taken 
on the James Clerk Maxwell Telescope in August 1997 
and reduced them using the SURF package. 

2.2. Infrared Data 

Cep A was observed on the University of Hawai'i 2.2m 
telescope at the f/10 focus with QUIRC (QUick InfraRed 
Camera) in K band (2.200 /im) in July and August 2003. 
The plate scale of the telescope is 0.1886 arcsec pixel -1 
and the field-of-view 193" x 193". The total on source 
integration time was 23.2 minutes. The average seeing was 



sion of the Metropolitan District Commission, Commonwealth 
of Massachusetts. 

3 The HHT is operated by the Submillimeter Telescope Ob- 
servatory on behalf of Steward Observatory and the Max- 
Planck-Institut fuer Radioastronomie. 
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Fig. 1. Grayscale K-band image over-plotted with SCUBA 850 /im contours of Cep A showing the stellar and dust 
content. The central position is at a(2000) = 22 h 56 m 18!9, 6(2000) = 62°01'42'.'6. The 850 contours start at 1 Jy 
beam -1 and increment by 2 Jy beam -1 . The large and small white plus signs show the location of the IRAS and MSX 
point sources respectively (each with a positional accuracy of less than 5"). The small white circle near the center 
of the map marks the location of the high mass protostar, HW2. The arrows show the direction of the bipolar CO 
outflow observed by |Rodriguez et al~| (1980). The dashed circle at the western edge of the map shows the location of 
the starless core seen in the CS maps and the dotted circle, centered at (10,10) with diameter 2' defines the region 
over which the spectral line averages were calculated. 



l'.'l FWHM. The observations were carried out by taking 
alternate object and sky exposures. All frames were flat- 
ficldcd using a normalized incandescent light dome flat. 
Sky frames were obtained by computing the median of 
several sky exposures close in time to a given object frame. 
The sky frames were subtracted from the object frames. 
The individual images were registered and co-added to 



produce a single image. Unfortunately, the conditions were 
non-photometric but we conservatively estimate that we 
should be able to detect to a limiting magnitude of at least 
19, which is 3.4 magnitudes fainter than 2MASS. 
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Fig. 2. Comparison of line and dust emission. The line maps have a resolution of 10" and are contoured with starting 
levels and increments equal to 4.0 K km s -1 for CS (2-1) and 5.0 for N2H 4 " (1-0). The dust emission has a resolution 
of 14" and is shown as a grayscale ranging from 1 to 15 Jy beam -1 . Symbols are as in Figure ^ 



3. Analysis 

3.1. Continuum and Line Maps 

Figure ^ shows contours of the 850 /im continuum emis- 
sion overlaid on the grayscale K-band image. Note that for 
35 K dust and the SCUBA beamsize of 14", 1 Jy beam -1 
corresponds to a visual extinction, Ay = 24, or an extinc- 
tion at K-band, Ak = 2.4. The infrared nebulosity lies 
within the extended sub-millimeter emission (Ak ^5 — 8) 
but is offset from the peak. There is one IRAS and one 
MSX point source in the region, slightly offset from each 
other. The IRAS-HiRes 12 /im and MSX 10 /im emis- 
sion are both associated with the near-infrared nebulos- 
ity, but the 60 /im IRAS-HiRes image peak is closer to 
the maximum of the 850 /im map, which traces the cooler 
dust. This suggests a spread in ages and location of star 
formation in Cep A, with the youngest protostars more 
deeply embedded in the core and invisible at near- and 
mid-infrared wavelengths. 

Spectral line maps and their comparison with the dust 
emission are shown in Figures [21 El In each case the emis- 
sion has been integrated over the full extent of the line, 
including all hyperfine components in the case of N2H + . 
Figure El shows the combined BIMA+FCRAO data at 
10" resolution. Figure shows these maps and the higher 
transition data smoothed to a uniform resolution of 30". 
Due to their high dipole moments, CS and N2H + both 
trace high volume densities over a range, riH 2 ~ 10 4 — 10 6 
cm~ 3 , for these transitions. The complex structure appar- 
ent in the line maps contrasts with the relative simplicity 



of the dust continuum emission. The spatial distribution 
of N2H" 1 " is similar to that of the ammonia UTorrelles et al.l 
1993) and generally follows the 850 /im continuum emis- 
sion, but N2H" 1 " is notably absent toward the core center 
where the CS is strongest. This behavior is opposite to the 
situation in low mass cores where N2H + is more centrally 
concentrated than CS (JTafalla et al.l 20021. 

The maps also reveal an apparently starless core lo- 
cated ~ 90" to the west of the bright 850 /im peak, indi- 
cated by a cross on Figure^ The core is detected in CS(2- 
1) and (5-4) but has only very weak N2H + (l-0) emission 
(Figure^ and no stars are seen in the near-infrared image. 
This core is on the edge of the SCUBA map where there 
is extended emission but no significant peak. The pres- 
ence of CS (5-4) emission indicates large enough densities 
that we would expect to observe N2H" 1 " emission, as in the 
large core. Since this emission is very weak, there must be 
a variation of abundances between the two cores. Poten- 
tial reasons for the wide range of abundances throughout 
the region are discussed in 21 

3.2. Radiative Transfer Modeling 

Our understanding of the structure and dynamics in Cep 
A is complicated by the small scale chemical variations 
within the core. Henceforth, we restrict attention to the 
large scale properties of the region as if observed with a 
2' gaussian beam centered at offset (10,10). The size and 
location of the averaging region was chosen so as to be 
broadly centered on the sub-millimeter continuum map 
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the red-shifted one, indicating that infall is the dominant 
effect on the scale of the core ( |Leung fc Brown| 1977). 

In order to examine the properties of the large scale 
structure and dynamics, we modeled the core average 
line profiles using the radiative transfer code, ratran 
pfogerheijde fe van der Tak | 2000). C ollisional rate coeffi- 
cient s for CS ar e from lTurner et all (1992 s ) and for N 2 H+ 
from Monteiro (1985 and references therein). From the 
brightness temperature of the N2H + and assuming an ex- 
citation temperature equal to the dust temperature, 35 K, 
we estimate the optical depth of each hyperfine compo- 
nent to be less than 0.1. For the purposes of the modeling, 
therefore, we can ignore the complicated hyperfine struc- 
ture (7 components at J — 1-0 and 45 at J = 3-2) and 
simply match the integrated intensity of each rotational 
transition. 

The inputs to the radiative transfer model include the 
density, temperature, and velocity structure of the core 
and the abundances of each molecule. However, many pa- 
rameters are constrained by related observations and re- 
sults. The remaining free parameters, eight in all, were 
then varied so as to fit the profiles of each of the three CS 
transitions and the integrated intensity of the two N2H" 1 " 
transitions. 



3.2.1. Model inputs 



Fig. 3. Line maps of the Cepheus A core shown as 
grayscale with contours overlayed. Positions are given in 
offset coordinates with the same central position as Fig- 
ure H CS maps are shown in the left panels and N2H+ 
in the right. The maps have been smoothed to a uniform 
30" resolution. The contour starting levels and increments 
for each map are 4.0 K km s" 1 (CS 2-1), 4.0 (CS 5-4), 
1.5 (CS 7-6), 5.0 (N 2 H+ 1-0), and 3.0 (N 2 H+ 3-2). Note 
the isolated core in the CS(2-1) map ~ 90" west of center 
(see text for discussion) . The 2' diameter circle shows the 
region over which the data were averaged to analyze the 
large scale dynamics of the core, and symbols are as in 
Figured 

but also to include the most prominent CS and N2H" 1 " 
structures within the core. The averaging region is indi- 
cated on the continuum and line maps in Figures ^ and 
[3J Averaged line profiles are displayed in Figure This 
figure also includes the line profile for C 34 S(2-1). This iso- 
topomer of CS, which is less abundant and therefore op- 
tically thinner, peaks near the same velocity as the dip in 
the average CS(2-1) spectrum. We therefore conclude that 
the two peaks in the latter are due to radiative effects (self- 
absorption) and not two features along the line of sight. 
Although individual CS profiles show both blue- and red- 
shifted self-absorption that may be due to a mix of infall, 
outflow, and rotation IjDi Francesco et al.1 20011. the aver- 
age CS(2-1) profile has a blue-shifted peak brighter than 



The 850 fim map in Figure^allows us to measure the col- 
umn density profile of the core and thereby estimate the 
volume density. Since the radiative transfer model is one- 
dimensional, we approximated the density structure as a 
radial function by calculating the average flux in concen- 
tric elliptical annuli centered on the peak of the emission. 
The equivalent radii were defined as the geometric mean 
of the semi-major and semi-minor axes of each ellipse. 
The column density was determined by assuming a dust 
temperature of 35 K ( |Moriarty-Schieve n et al.| 1991) and 
a mass-opacity n = 0.02 cm 2 g~ L ( |Ossenkopf fe Henning| 
1994). The volume density was then derived by assuming a 
path length through the core equal to twice the equivalent 
radius. The resulting density profile was fit to a Plummcr- 
like model, 



n H 2 ( r ) = "0 



-a/2 



where n = 1.1 xlO 7 cm~ 3 , r = 0.02 pc and a — 2.0. The 
model discretizes this density profile over nine logarithmi- 
cally spaced shells. Both measured and discretized density 
profiles are displayed in Fig. there are no data points 
for radii smaller than 0.03 pc, the equivalent radius corre- 
sponding to a semi-minor axis of half the 850 /im beam- 
size, and for radii larger than 0.18 pc, which corresponds 
to the extent of the 850 fim map. 

Additional inputs to the model include a tempera- 
ture of 35 K, derived from a graybody fit to the SED 
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Fig. 4. Averages of CS and C 34 S (2-1) (left) and N 2 H+ 
(right) spectra over the central core region outlined by 
the circles in Figures 1 and|3| The heavy solid lines show 
the observed data and the dashed line shows the modeled 
profiles for the CS lines. The multiple peaks in the N 2 H + 
spectra is due to hyperfine structure which was not mod- 
eled (the fits matched only the integrated intensity of each 
rotational transition) . 

( |Moriarty-Schieven et al.| 1991), and a systemic velocity 
and velocity dispersion derived from fitting the optically 
thin N 2 H + (l-0) line, v COIC = —10.5 km s _1 and a COIC = 
1.2 km s _1 , respectively. A constant temperature and ve- 
locity dispersion were adequate to model the core. 

To model the observed self- absorption in the CS(2- 
1) line, we found that we required a low excitation (and 
therefore low density) and low velocity dispersion (<r env = 
1.1 km s _1 ) outer shell. The low velocity dispersion is re- 
quired by the narrowness of the observed absorption dip. 
The overall model has 10 shells therefore, but it can ef- 
fectively be considered as two layers; a power law core 
on the inside and a low density outer envelope. The free 
parameters in the model are the inner core cutoff radius, 
the size, density, velocity dispersion and relative velocity 
of the envelope and the molecular abundances in the core 
and envelope. 

The values for the molecular abundances were guided 
by observations of other cores and theoretical mod- 



els. In cold molecular cloud cores prior to star forma- 
tion, the chemistry is dominated by low-temperature 
gas-phase ion-molecule and neutral-neutral reactions 
Ijvan Dishoeck fc Blakel 1998). During the cold collapse 
phase, however, the density becomes so high that many 
molecules freeze onto grain surfaces. ITafalla et al.l (1998) 
find that, in low mass star forming cores, the abundances 
of the tightly bound sulfur-bearing molecules such as CS 
begin to exhibit large depletions at densities in the range, 
tih 2 ~ 2 - 6 x 10 4 cm~ 3 . This behavior is in contrast to 
that of N 2 H + which, due to the low binding energy of the 
precursor molecule N 2 , depletes only at the highest den- 
sities (Bcrgin & Lange^l997). Based on these results, we 
assume a constant abundance of N 2 H + but allow the CS 
abundance to vary. As we show later, we fit the data with 
a simple "jump" model where the CS is depleted in the 
inner dense core relative to a lower density outer envelope. 
This is a similar abundance profile to the detailed chemi- 
cal modeling of the high mass star forming region, AFGL 
2591, by |Doty et at] (2002). 

3.2.2. Model results 

By varying the density, size and velocity of the envelope, 
and the CS abundance in the core and envelope, we were 
able to fit the integrated intensities of the three observed 
CS transitions to within 20% and reproduce the asymme- 
try in the (2-1) line reasonably well. The observed and 
model CS spectra are compared in Figure [3J The param- 
eters of the fit are tabulated in Tabled where r coro is the 
inner core cut-off radius, r onv , n cnv and v- ln — v cnv — u core 
are the size, density and relative velocity of the outer enve- 
lope, and x(CS) corc and a;(CS) cn v are the CS abundances 
in the core and the envelope respectively. The density, 
velocity, and CS abundance profiles are graphed in Fig- 
ure and annotated with the above parameters. Simulta- 
neously, the N 2 H + (1-0) and (3-2) integrated intensities 
were also matched: we found 8.1 and 9.5 K km s _1 re- 
spectively for the modeled values, whereas the observed 
values are 9.7 and 8.2 K km s _1 for the (1-0) and (3- 
2) transitions respectively. These results correspond to a 
match of 16 and 15% respectively. Finally, we also fitted 
the C 34 S (2-1) integrated intensity, using an abundance 
a;(C 34 S) = z(CS)/30. We found a modeled value of 1.4 K 
km s _1 which is within 8% of the observed value (1.5 K 
km s _1 ). 

Despite the good overall fit, some slight discrepancies 
remain. The observed CS(5-4) is stronger than the model 
by 20%. This may be due to calibration error or may re- 
flect a more complex CS abundance profile in the core. 
The observed CS(2-1) profile is slightly narrower than the 
model spectrum, perhaps due to the constant linewidth 
assumption in the core (the velocity dispersion of cores is 
expected to decrease with radius. IL"arsonl 1981). A more re- 
fined model with additional parameters would fit the data 
more closely but probably not with greater significance. 
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Table 1. Fit parameters. 
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Fig. 5. Radial profiles of H2 volume density, core-envelope 
CS abundance and velocity profiles used in the best fit 
model. Relative motions in the inner region, r < r core , 
were not well constrained by the data. The thin solid line 
shows the H2 density profile derived from the 850 fim data 
(see text). 

Since the model reproduces the blue-red peak asym- 
metry in the CS(2-1) spectrum well we are confident that 
we have accurately measured the average infall speed of 
the envelope, V[ n — —0.22 km s" 1 . However, the size of 
the infalling region is not well determined because the av- 
erage profiles of the higher transition CS lines are not self- 
absorbed and therefore relative motions of the higher den- 
sity gas could not be constrained. Nevertheless, we were 
unable to model the spectra with a static outer shell and 
an inner collapse: no such set of parameters could repro- 
duce the dip in the CS (2-1) spectrum. That is, our results 
indicate a large scale collapse from the outside-in. 

4. Discussion and Summary 

The line observations reveal a complex chemistry in the 
core. The CS emission is concentrated toward the cen- 
ter near the peak of the 850 (im dust emission and the 
youngest, most embedded protostars. However, the N2H + 
map shows two prominent cores offset on either side of the 
dust peak. The presence of CS and absence of N2H + to- 
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ward the star forming center of the core is likely due to the 
fact that neutral molecules released from the dust grains 
in the hotter region surrounding the protostars preferen- 
tially destroy ions such as N2H 4 " ( |Ber gin 2000). 

The maps also reveal a small starless core toward 
the west of the main core. The relatively strong CS 
and weak N2H+ emission toward this core suggests that 
it has only recently formed (Bcrgin & Langcr 1997). 
Wil liams fc Myers| (1999) found a starless core with simi- 
lar chemical properties in the Serpens NW cluster 

Despite the chemical complexities on small scales, the 
average CS(2-1) spectrum is asymmetrically self-absorbed 
suggesting large scale collapse. Using the Hogerheijde & 
van der Tak radiative transfer code, we fit the average 
CS and N 2 H+ spectra with a spherical model consist- 
ing of an inner region with a Plummcr-like density pro- 
file measured from archival SCUBA 850 /im data with 
constant temperature and linewidth. The best fit model 
that matches the three CS profiles has a low CS abun- 
dance in the inner region and an outer, infalling enve- 
lope with a low density and higher CS abundance. The 
depletion toward the center matches chemical evolution 
model expectations | |Bergin fe Langer| 1997) and the en- 
velope CS abundance is similar to that in the extended 
ridge of Orion Ijvan Dishoeck fc Blakel 1998). The fit also 
matched the integrated intensities of the two N2H" 1 " spec- 
tra with a constant abundance similar to that found in 
B68 ( |Bergin et al.| 2002). In practice the maps show that 
the N2H + must deplete toward the core center (see also 
|Doty et al.| 2002) but we have not attempted to match 
the complicated small scale chemistry and our model fit 
for the N2H + abundance should only be considered an 
average, weighted by column density, over the core. 

The CS(2-1) self-absorption requires a large scale 
outside-in collapse. The velocity of the collapse could be 
accurately measured but the depth of the collapse region 
could not, due to the absence of self-absorption in the 
higher transition CS lines. At 0.22 km s" 1 , the infall ve- 
locity is sub-sonic for a gas temperature of 35 K. The mass 
infall rate of the envelope can be estimated from the ratio 
of its mass, M env = 240 M Q determined from its size and 
density, and the time for the outer edge to reach the core 
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center. 

ii V -^env^in — 1 r\ — 5 ti /r — 1 

M in = = 7 x 10 M yr \ 

^env ~r 7"corc 

This may only be a lower limit to the total mass infall 
rate at the center if the inner region is also collapsing. 
Nevertheless, the envelope mass infall rate alone is more 
than an order of magnitude higher than typical mass infall 
rates for solar mass protostars l|Zhoul l995h 

Our data do not rule out inside-out collapse motions 
around individual protostars at higher densities on smaller 
size scales since the r = 1 surface of the CS(2-1) emission 
occurs at low densities and therefore at large scales. The 
multitude of sources, powerful outflows, and the complex 
chemistry would likely make an investigation of the small 
scale motions around individual protostars quite challeng- 
ing. 

On large scales, however, our picture is of a core with 
a Plummer-like density profile accreting low density gas 
sub-sonically. The CS abundance in the infalling envelope 
is similar to undepleted values in the ISM. The dissi- 
pation of turbulent support resulting in "cooling flows" 
may lead to core growth in this manner iTSTakanpl 1998 
; |Myers fc Lazarian| 1998 ; |WflI iams &: Myers| 2000). The 
growing availability of dust continuum maps and multi- 
transition, multi-species line observations will lead to more 
refined structural and dynamical modeling and compar- 
isons between different star forming environments in the 
future. 
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